p=[1 9 135/4 135/2 1215/16 729/16 729/64];
x1=roots(p)
p1=poly2sym(p);
x2=solve(p1)